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ABSTRACT 

Context. It is well established that between 380000 and 1 billion years after the Big Bang the Inter Galactic Medium 
(IGM) underwent a "phase transformation" from cold and fully neutral to warm (« 10 4 K) and ionized. Whether this 
phase transformation was fully driven and completed by photoionization by young hot stars is a question of topical 
interest in cosmology 

Aims. We propose here that besides the ultraviolet radiation from massive stars, feedback from accreting black holes in 
high-mass X-ray binaries (BH-HMXBs) was an additional, important source of heating and reionization of the IGM in 
regions of low gas density at large distances from star-forming galaxies. 

Methods. We use current theoretical models on the formation and evolution of primitive massive stars of low metallicity, 
and the observations of compact stellar remnants in the near and distant universe, to infer that a significant fraction 
of the first generations of massive stars end up as BH-HMXBs. 

Results. The total number of energetic ionizing photons from an accreting stellar black hole in an HMXB is comparable 
to the total number of ionizing photons of its progenitor star. However, the X-ray photons emitted by the accreting black 
hole are capable of producing several secondary ionizations and the ionizing power of the resulting black hole could be 
greater than that of its progenitor. Feedback by the large populations of BH-HMXBs heats the IGM to temperatures 
of ~ 10 4 K and maintains it ionized on large distance scales. 

Conclusions. BH-HMXBs determine the early thermal history of the universe and mantain it as ionized over large 
volumes of space in regions of low density. This has a direct impact on the properties of the faintest galaxies at high 
redshifts, the smallest dwarf galaxies in the local universe, and on the existing and future surveys at radio wavelengths 
of atomic hydrogen in the early universe. 
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1. Introduction 

The so-called "dark ages" of the universe started « 380000 
years after the Big Bang as matter cooled down and space 
became filled with neutral hydrogen. This phase of the uni- 
verse lasted up to about a billion years, when the com- 
plex process of reionization of the IGM was completed 
(jLoeb 2010[) . It is currently believed that most of the reion- 
ization was caused by the ultraviolet radiation from massive 
stars formed in the first generations of galaxies. However, 
it is uncertain what fraction of ionizing ultraviolet pho- 
tons could escape from primitive galaxies to produce and 
maintain the ionization far from galaxies in low density 
regions of the IGM. Recent observations with the Hubble 
Space Telescope suggest that the rest-frame ultraviolet ra- 
diation from the most distant galaxies detected so far at 
the heart of the dark ages is not enough to heat and ionize 
the IGM over large volumes o f space. To solve this appar- 
ent "photon-starved problem" (Bouwens et al. 2010|) . it has 



been suggested by Lehnert et al. (2010) that a fainter pop- 
ulation of galaxies below the present detection limit could 
contribute significantly to the reionization. 

X-rays from accreting black holes have a longer mean 
free path than the ultraviolet photons from massive stars. 
In this context Madau et al. (2004) and Ricotti & Ostriker 
(2004) have suggested that a smaller fraction of ioniz- 
ing photons were provided by primordial black holes of 
intermediate mass ( "miniquasars" , at z > 10) accret- 
ing by Bondi-Hoyle from the surrounding gas. However, 
feedback from Bondi-Hoyle accretion to solitary black 
holes significantly suppresses both any further inflow 
dAlvarez, Wise k Abel 2"009{|Milosavljevic et al. 2009[ ) and 
the consequent injection of radiation and high-energy par- 
ticles to the surrounding medium. 

We propose that BH-HMXBs at z > 6, namely, 
the fossils of massive stars are an important -so far 
overlooked- agent in the complex process of the reion- 
ization of the universe. In the context of the cur- 
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rent models on the formation (Krumholz et al. 2009: 
|Turk, Abel fc O'Shea 2009| |Stacy, Greif fc Bromm 2010P 



and collapse QHeger et al. 2003[ |Meynet fc Maeder 2005 
|Georgy et al. 2009[ ILinden et al. 2010|) of primordial stars 
an attractive and realistic alternative to the hypothesis of 
quasi-radial, Bondi-like accretion on solitary black holes of 
intermediate mass ( "miniquasars" ) is accretion on stellar 
black holes from high-mass stars in binary systems, namely, 
"microquasars" ( Mirabel fc Rodriguez 199 9). As shown be- 
low, the formation rate of BH-HMXBs must have been very 
large in the young Universe, playing an important role in 
the thermal history of the IGM, and a complementary role 
to that of their progenitor stars in the re-ionization process 
of the IGM over large volumes of space. It is this scenario 
that we investigate here. 



2. Cosmic evolution of BH-HMXBs: A prediction 
from current theoretical models 

Recent hydrodynamic simulations of the formation of the 
first generations of stars show that a substantial fraction 
of stars in primordial galaxies form as binaries with typi- 
cal masses of tens of solar masses (Krumholz et al. 2009; 
Turk, Abel fc O'Shea 20091 |Stacy, Greif fc Bromm 20T0l ). 



Models of single stars with very low metal content and ini- 
tial masses of a few tens of solar masses show that they 
collapse directl y with no energetic natal kicks, and end 
as black holes flHeger et al. 2003[ |Meynet fc Maeder 2005] 
|Georgy et al. 2009| ). 

On the other hand, a recent model of binary evolu- 
tion of massive stars by Linden et al. (2010) show that 
the number of HMXBs and ULXs, their time evolution, 
and orbital period distribution, are strongly metallicity de- 
pendent. Linden et al. (2010) find that ULXs formed in a 
typical starburst of 10 6 Mq with Z = 0.02 Zq outnumber 
ULXs formed with Z = Zq by a factor of 5, and after 10 
Myr by almost three orders of magnitude. Besides, at Z = 
0.02 Zq, among the ULX population, > 95% of the com- 
pact objects are black holes formed by direct collapse and 
therefore after black hole formation it remains gravitation- 
ally bound to a companion, donor star. Most of the orbital 
periods at Z = 0.02 Zq are less than 3 days and accre- 
tion is by Roche lobe overflow which creates very luminous 
and persistent BH-HMXBs. Probably, this trend continues 
for starbursts with the metallicities (Z < 0.02Zq) of the 
reionization era. 

These models imply that the majority of the first gener- 
ations of high mass stellar binaries remain gravitationally 
bound after the formation of black holes. Massive stellar 
binaries can thus become BH-HMXB microquasars, which 
are sources of UV photons, X-rays, massive winds, and rela- 
tivistic jets (Mirabel & Rodriguez 1999). Therefore, in the 
context of the models of massive stellar evolution and the 
cosmic evolution of metallicity it is expected that: 1) the 
fraction of black holes to neutron stars and 2) the fraction 
of black hole binaries to solitary black holes, should increase 
with redshift. That is, the rate of formation of bright BH- 
HMXBs was likely much larger in the early Universe than 
at present. 



3. Formation rate of stellar black holes as a 
function of metallicity: Observations 

The cosmic evolution of BH-HMXBs inferred from theoret- 
ical models is consistent with the following observational 
studies of stellar black holes and neutron stars in the near 
and distant universe: 

1. The mass of black holes in HMXBs seems to be 
a decreasing function of the host galaxy metallicity 
( Crowthcr et al. 2010l and references therein). The black 
holes in the binaries M 33 X-7, NGC 300 X-l, and IC10 
X-l are in low-metallicity galaxies and have masses - 
determined dynamically- in the range of 16 to 30 so- 
lar masses, which are higher than the mass of any 
known stellar compact source in the Milky Way and 
Andromeda galaxies, which have higher metallicities. 
However, while the model by Linden et al. (2010) 
supports the formation of HMXBs and ULXs in low- 
metallicity environments, they conclude that it is diffi- 
cult to create very massive black holes through common 
envelope phases, since this tends to strip a high fraction 
of the primary envelope. Given the low number statistics 
of the known dynamic masses of black holes in HMXBs, 
it is possible that the relatively high masses of 16 to 30 
solar masses come from the selection of the brightest 
sources, namely, those that are at the tip of the iceberg. 

2. It is believed that the majority of ultraluminous X-ray 
sources (ULXs) found in external galaxies are HMXBs 
that contain black holes accreting at Super-Eddington 
rates ( [Gladstone, Roberts fc Done 2 009). In fact, the 
occurrence rate per unit galaxy mass of ULXs observed 
in nearby galaxies is a decreasing function of the galaxy 
mass - hence of the metallicity - of the host galaxy 
dZamperi fc Roberts 2009[ ). 

3. The space kinematics of Galactic X-ray binaries that 
contain black holes with more than ten solar masses 
provides evidence of black hole formation by implo- 
sion, with no large kicks due to energetic supernovae 
( |Mirabel fc Rodrigues 2003} IMirabel 2010|) . 

4. Observations now support the notion that massive 
stars with high metal content may end as neutron 
stars instead of bl ack holes (|Meynet fc Maeder 2005 



IGeorgy et al. 2009[ ILinden et al. 2010p . 



Recently 

formed neutron stars observed as soft gamma ray 
repeaters and anomalous X-ray pulsars are found in 
young clusters of large metal content that contain stars 
with masses of 40-50 solar masses QFiger et al. 2005] 
IMuno et al. 2006|) . 

5. It is believed that the majority of gamma ray bursts 
of long duration (LGRBs) mark the formation of black 
holes by the collapse of massive stars. Although a 
fraction of dark LGRBs may require local extinction 
columns of Av > 1 mag, the majority of their hosts 
are faint, irregular galaxies with global limited chemical 
evolution (|Le Floc'h et al. 20031 IFruchter et al. 2006| 
IHan et al. 20T0l |Levesque et al. 2010fl . The proper- 
ties of GRB 090423 at z = 8.1 are similar to 
those of GRBs observed at low/intermediate redshifts 
(|Salvaterra et al. 20 09). suggesting that the mecha- 
nisms and progenitors that gave rise to this burst are 
not too different from those producing GRBs with iden- 
tified hosts. 

6. There is increasing evidence for an enhanced LGRB 
rate at z > 3 (Daignc, Rossi & Mochkovitch 2006; 
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Qin 2010), as expected from the increase in the 
specific sta r formation rate (S FR) with decreasing 
metallicity (Man nucci et al. 2010|) . 



4. Ionizing power of a stellar black hole in an 
HMXB relative to its progenitor star 

In the following we compute the number of ionizing soft 
X-rays and UV photons from the accretion disk of a BH- 
HMXB, and then compare its ionization power with that of 
its progenitor massive star. To this end we assume that the 
black hole of mass Mbh is accreting at a fraction f e dd of its 
Eddington luminosity for a time t acc . Then the ratio of the 
total number of ionizing photons emitted by the accreting 
black hole to that emitted by the progenitor star is given 
by 
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where N p h ot denotes the number of ionizing photons emit- 
ted per hydrogen nucleus involved in star formation, (Ej) 
denotes the mean photon energy of the radiation emitted 
by the accreting black hole, and f CSC: * (f eS c,BH) denotes the 
fraction of ionizing photons emitted by the star (accreting 
black hole) that escape from the galaxy and contribute to 
the heating and reionization of the IGM. 

We have substituted reasonable numbers for each of the 
parameters: N p h o t=64000 corresponds to metal-free form- 
ing stars with a top-heavy Initial Mass Function (IMF). 
For a normal stellar population this number would be 16 
times lower (Schacrer 2003}. The escape fraction of ion- 
izing photons emitted by stars was taken to be f eS c,* 
0.1. which is consistent with the mean observed escape 
fraction of ionizing photons from star-forming galaxies 
( Shapley et al. 2006| ) at z=2-3. 

HMXBs are expected to inject photons into the IGM at 
a rate close to Eddington (and possibly super Eddington), 
and our choice f c dd — 0.1 is likely to be conservative 
(see footnote 1 ). The accretion lasts t acc w 20 Myrs, 
a mean lifetime of donor stars of M» = 10 — 30M(T) 
dTurk, Abel fc O'Shea 20091 IStacy, Greif fe Bromm 2010[ ). 
The mean photon energy of the radiation emitted by the ac- 
creting source depends on the assumed spectral shape, but 
(E y ) — 1 keV appears -within a factor of a few- of values 
derived for both the thermal and power-law components 
of the black hole spectra. Finally, f eS c,BH = 1.0 is not well 
constrained. We took f esc ,BH to be larger than f esCi * simply 
because energetic photons can propagate much more easily 
through HI column densities in the range 10 17 — 10 20 cm~ 2 . 
For reasonable choices for each of these model parameters, 
we find that an accreting black hole in a high-mass binary 
emits a total number of ionizing photons that is comparable 
to its progenitor star. 

However, it should be kept in mind that the ioniz- 
ing photons emitted by the accreting black hole are more 
energetic than those from the progenitor star and ca- 
pable of ionizing more than one hydrogen atom. In a 
fully neutral medium, the number of secondary ioniza- 
tions N see = 25(E 7 /lkeV), wher e E 1 is the photon en- 
ergy flShull fc van S tecnbe rg 1985] ). Therefore, the ionizing 



power of the resulting black hole could be greater than that 
of its progenitor. 



5. BH-HMXBs and massive star formation rates in 
the epoch of reionization 

The progenitor star of the black hole should have formed 
in a molecular cloud, which in turn may have formed more 
stars, each of which contributed to the total number of ion- 
izing photons emitted by stars. Therefore, in the following 
we look at the emission of ionizing radiation from a star- 
forming region as a whole. 

Observations of galaxies in the local universe show that 
their X-ray luminosity in the energy range 2-10 keV cor- 
relates strongly with the rate at which they are forming 
stars ( |Grim, Gilfanov Ik Sunyaev 2003[ ) . This local correla- 
tion states that the X-ray luminosity (in erg s _1 ) scales 
with the SFR in M yr" 1 , as L 2 - W = 7 1 39 SFR erg s _1 . 
When modeling the impact of X-ray emission from galaxies 
at very high redshift, the following more general correlation 
is used (jFurlanetto 2006|) : 



L 2 -io = fx x 3.5 1 40 x SFR erg s" 1 . 



(2) 



Here, the parameter fx accounts for the likely case 
that the normalization of the observed correlation depends 
on redshift. Observations indicate that fx = 0.2 for local 
galaxies. This important parameter fx depends on several 
physical processes (all of which are expected to change with 
redshift). To illustrate this quantitatively, we express the 
X-ray luminosity of a star-forming galaxy in more funda- 
mental quantities as 



£2-10 — 

12-10 X X tacc X f bin X f] 

or 



? E dd x 1.510 38 erg s^M" 1 ( 3 ) 



L2-10 — 12-10 x ibh x SFR 



x t a 



IBH a or iv /n u acc , , 

xfbin x f Edd x 1.5 1 38 erg s^M" 1 W 

Most of the parameters in this equation were intro- 
duced earlier. The fraction 12-10 denotes the fraction of 
the total luminosity that emerges in the 2-10 keV band. 
For the power law component of the spectra we expect 
12-10 = 0.1 — 0.5 (equal flux per logarithmic bin of en- 
ergy or steeper), and we conservatively adopt 12-10 = 0.1. 
We also introduced the parameter fsH, which relates black 
hole formation rate dMs^/dt to SFR. This fraction can 
be computed for a given initial IMF under the assump- 
tion that every star above some critical mass M*. cr ;t ends 
up as black hole. As argued previously, due to mass lost 
by metallicity-dependent stellar winds even massive stars 
may end up as neutron stars instead of directly as black 
holes. However, in the evolution of a close massive binary 
of low metallicity, due to mass transfer and common enve- 
lope phase the primary could lose mass and end its life as 
a neutron star rather than a black hole, leading to suppres- 
sion of black hole formation. Linden et al. (private commu- 
nication) discussed the importance and frequency of such 
scenario finding that the transition from neutron star to 
black hole dominated HMXBs is very sharp and metallic- 
ity dependent, and that for Z < 0.02Zq all primary stars 
with M > 2OM in HMBs end as black holes in HMXBs. 
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Therefore we assume here that mass transfer in close bi- 
naries has little impact on the mass range of stars becom- 
ing black holes. For the conservative choices of a Salpeter 
IMF in the range M low = Q.IMq, M up = 100M Q , and 
M» jCr it = 25Mq we find {eh = 0.03. Finally, we expect 
black holes in binaries to be very efficient accreters, while 
isolated black holes are not. The parameter f^ denotes the 
mass fraction of accreting black holes in binaries. The total 
mass fraction of massive stars in binaries is close to unity, 
and for simplicity we assume that this total mass fraction is 
50 %. Since BH-HMXBs are persistent sources of radiation 
and jets for our model we assume %in = 0.5. 

Both Eqs. (2) and (4) depend linearly on SFR, and we 
can express the parameter fx as a combination of physical 
parameters: 




Linden ct al. (2010) estimate that at ULX luminosities 
the number of Z = 0.02 Zq sources should outnumber the 
Z = Zq HMXBs by a factor of 5. Therefore, in our fiducial 
conservative model for primordial starbursts fx is at least 
one order of magnitude higher than the locally observed 
value flGrim, Gilfanov k Sunyaev 2003[ ). There are several 
reasons why we expect fx to be much higher in the young 
Universe. (1) Microquasars can have spectra that are harder 
and 12--10 could be in some cases as high as 12-10 = 0.5; and 
(2) ign is higher in low metallicity environments, i.e., the 
formation of black holes in metal-enriched environments of 
galaxies in the local universe is likely to be strongly sup- 
pressed compared to more pristine environments. (3) Stars 
in the young Universe wre very likely formed with an IMF 
that was more top-heavy. This evolution in the IMF alone 
could boost {bh to much higher values than the fiducial 
assumed value in Eq. (4). Finaly, the IGM at z = 10 would 
essentially be transparent to the hard X-ray photons (> 1 
keV), and following Dijkstra, Haiman k Loeb (2004) we 
explicitly verified that these models are consistent with the 
observed unresolved soft X-ray background. 

1 To estimate the ionizing power of BH-HMXBs in the early 
universe we use the observed spectra of the Galactic black hole 
binary Cygnus X-l, which has a blue supergiant companion, as 
template for moderate accretion rates. In this system, accretion 
is persistent and composed mainly by the donor stellar wind. 
Its X-ray spectrum is characterized by two components: a UV- 
soft X-ray bump due to thermal emission from the accretion 
disk with a typical temperature of as 7 eV, and a non-thermal 
power-law component of hard X-rays that result from Compton 
up scattering of thermal accretion disk and synchrotron pho- 
tons by a hot coronal plasma and/or jet. In addition to X-rays, 
these sources can indeed produce jets and winds of accelerated 
particles, which in turn may heat and ionize the surrounding 
medium. 

On the other hand, for higher accretion rates one can use the 
ULXs observed in the local universe as templates of BH-HMXBs 
in the early universe. ULXs often have spectra that resemble the 
very high (super-soft) state observed in some Galactic black hole 
binaries (e.g. GRS 1915+105). The majority of ULXs exhibit a 
complex curvature, which can be modeled by a cool disk compo- 



6. Stellar black holes and the thermal history of 
the IGM 

It is an open question how significant may be boosting 
the X-ray emissivity of star-forming galaxies in the high- 
redshift Universe for the global ionized fraction (i.e. aver- 
aged over the entire volume of the observable Universe). 
However, it has been shown ([Furlanetto 2006) that soft 
X-rays (E 7 < 2 keV, IPritchard k Furlanetto 2007|) and 
inverse-Compton scattering from relativistic electrons could 
have profound implications for the amount of heating of the 
low-density neutral IGM. BH-HMXBs are powerful sources 
of soft X-rays and relativistic jets, and their heating could 
in turn affect the overall reionization process indirectly in 
ways that will need to be investigated further. In Fig [T] 
we show that fx (as defined in equations 5 and 6) is the 
parameter that determines the thermal history of the IGM. 
Increasing fx causes the neutral IGM to be heated earlier. 
As shown in the thermal evolution for fx = 10.0, when 
the gas temperature approaches 10 4 K, cooling through 
collisional excitation of the atomic hydrogen becomes ef- 
ficient, and further heating is not possible in practice. As 
discussed in section 8, the formation of low-mass galaxies 
in the neutral IGM at z=10-20 is suppressed when the gas 
temperature is as high as 10 4 K, namely, for fx > 5. For 
further details on this particular model we refer the reader 
to Pritchard k Furlanetto (2007). 

7. The 21cm line of HI during reionization 

The precise temperature evolution of the neutral IGM is 
known to strongly affect the global 21 cm signature ex- 
pected from neutral HI during the epoch of reionization 
dFurlanetto, Oh k Briggs 2 006). In Fig [2) we plot the evo- 
lution of the brightness temperature (averaged over the en- 
tire sky) as a function of redshift following the same pre- 
scriptions as in Pritchard k Loeb (2010), again for the same 
three models with different values of fx- The gas temper- 
ature couples to the excitation temperature through colli- 
sions, and through scattering of Lyman a photons. When 
this excitation temperature, also known as the spin tem- 
perature, is higher (lower) than the temperature of the 
cosmic microwave background (CMB), then it is possi- 
ble to observe hydrogen atoms in emission (absorption) 
against this CMB. The difference in spin and CMB tem- 
perature is referred to as the "brightness" temperature 
( jFurlanetto, Oh k Briggs 2006). Clearly, when the gas is 
heated earlier, the spin temperature can increase earlier, 
and hydrogen can be seen in emission earlier in the evolu- 
tion of the Universe. Increasing fx also reduces the interval 
in redshift - and therefore - frequency over which hydrogen 
can be seen in absorption against the CMB (corresponding 



nent, together with a power law that breaks above 3 keV, prob- 
ably due to a cool, opticall y thick corona produced by super - 
Eddington accretion flows ([Gladstone, Roberts & Done 2009). 
Examples of steady super-Eddington sources are also SS 433 in 
the Milky Way, which is blowing the nebula W50 laterally, and 
the microqu asar that is inflating the neb ula S26 in the galaxy 
NGC 7793 dPakull, Soria fc Motch 20To[ ). SS 433 injects more 
than 10 39 erg s~ in the interstellar medium and the microquasar 
in S26 more than 10 40 erg s _1 . The overall energy injected by 
these microquasars during their whole lifetime can be more than 
10 54 er(/, which is orders of magnitude more than the photonic 
and baryonic energy from a typical core collapse supernova. 
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Fig. 1. Thermal history of the low-density neutral inter 
galactic medium (IGM) due to heating by accreting stel- 
lar black holes in high mass X-ray binaries (BH-HMXBs) . 
This figure shows the gas temperature of the IGM as a func- 
tion of redshift z for three possible values of fx as defined 
in Eq. (6): fx = 0.1 [blue dashed line], fx = 1.0 [red dotted 
line], fx = 10.0 [black solid line]). As discussed in the text, 
most likely fx > 1- 



to negative brightness temperature). Single dipole exper- 
iments such as EDGES are currently attempting to mea- 
sure this signal ( |Bowman, Rogers fc Hewitt 2008[ ). One of 
the present challenges in observational astronomy is to di- 
rectly observe this 21 cm signal from neutral hydrogen in 
the young Universe ( |Morales fc Wyithe~2 010). This can be 
accomplishe d with single radio dipole experi ments, such 
as EDGES ( |Bowman, Rogers fc Hewitt 2008] ), which are 
potentially capable of d etecting the global 21 cm signal 
(jPritchard fc Loeb 2010|) at redshifts z < 30, directly mea- 
suring early heating of the IGM. 

A boosted X-ray emissivity of star forming also 
affects the fluctuations in the 21 cm background 
(jPritchard fc Furlanetto 20071 IPritchard fc Loeb 2008|) . 
Detecting these fluctuations is one of the prime scientific 
drivers for the next generations of radio interferometers 
such as the MWA, LOFAR, and SKA. The fluctuations in 
the 21 cm background radiation contain the largest amount 
of cosmological information ( |Loeb fc~Z aldarriaga 2004[ ), 
and are therefore invaluable constraints on cosmological 
parameters (as well as fundamental physics). However, "as- 
trophysics" introduces additional fluctuations, for example, 
through temperature fluctua tions and through fluctuations 
in the Lyman a backgro und (jPritchard fc Furlanetto 2007; 
IPritchard fc L oeb 2008). Both are sourced (on different 
scales) by galaxies that themselves provide biased tracers of 
the underlying density field. To fully exploit the rich data 
set provided by the 21 cm fluctuations therefore requires 
a good understanding of the nature of the astrophysical 
sources illuminating, heating, and ionizing the hydrogen 
in our Universe. Interestingly, one heats the neutral IGM 
earlier by boosting fx, which implies that temperature 
fluctuations are suppressed during the later stages of 
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Fig. 2. Brightness temperature of the hyperfine transition 
of the ground state of atomic hydrogen (the wavelength of 
this transition is 21 cm), averaged over the entire sky as a 
function of redshift z for the same three different values of 
fx as in Fig. [I] 



reionization, which improves the prospects for extracting 
cosmological information from the 21 cm background 
(jPritchard k. Furlanetto 20071 IPritchard fc Loeb 2008 j) . 



8. The role of stellar black holes in the formation 
of dwarf galaxies 

The cold dark matter model of the universe provides the 
framework for significant progress in understanding the 
large-scale properties and physical principles that govern 
the large-scale evolution of the universe during the first 
400 thousand years. However, it is still poorly understood 
how the first stars and black holes in galaxies were formed 
and the way that in less than a billon year these pristine 
objects re-ionized and re-heated most of the matter in the 
universe over large volumes of space. The apparent dispar- 
ity between the number of dwarf galaxies predicted by the 
cold dark matter model of the universe and the number of 
small galaxies observed so far in the halo of the Galaxy is a 
subject of topical interest in cosmology (jLoeb 2010|) . Power 
et al. (2009) had already pointed out the possible implica- 
tions of X-ray binaries in primordial globular clusters for 
the reionization, hence, galaxy formation at high redshifts. 

It is believed that the first stars had formed in gas clouds 
with a virial temperature of a few hundred K due to cooling 
by molecular hydrogen, H2 (jLoeb 2010 j) . But the UV radi- 
ation produced by these stars could have easily dissociated 
H2, making atomic hydrogen (H I) cooling necessary for fur- 
ther star formation. The galaxies that reionized the IGM 
were therefore likely to have a virial temperature above the 
H I cooling threshold of 10 4 K. X-ray, and UV heating by 
BH-HMXBs of the diffuse IGM during reionization would 
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have resulted in an additional increase in the total mini- 
mum galaxy mass M m i„ 

^n = 10 9 ( I5 ^)"'(^)- § 

(^) a (^)" § Af , 

where pc is the critical mass density for a flat Universe, 
p the mass density in the galaxy, p the mean molecular 
weight, z the redshift, and T the temperature of the IGM. 
Once the IGM was heated to a temperature of 10 4 K, dark 
matter halos with masses below 10 9 Mq could no longer 
accrete IGM material because the temperature of the in- 
falling gas increased by an extra order of magnitude as its 
density increased on its way into these galaxies. In that 
regime, only gaseous halos with virial temperatures above 
10 5 K could have accreted fresh IGM gas and converted it 
to stars. The census of dwarf galaxy satellites of the Milky 
Way requires a related suppression in the abundance of 
low-mass galaxies relative to low-mass dark matter halos 
(see IMunoz et al. 20091 and references therein) . The ther- 
mal history of the IGM therefore has a direct impact on 
the properties of the faintest galaxies at high redshifts, as 
well as on the smallest dwarf galaxies in the local universe. 

It is interesting to note that black holes of dif- 
ferent mass scales play a role in galaxy formation. 
Feedback from supermassive black holes halt star forma- 
tion, quenching the unlimited mass growth of massive 
galaxies (Cattaneo et al. 2009), and we show here that feed- 
back from stellar black holes in HMXBs during the reioniza- 
tion epoch suppress the number of low-mass dwarf galaxies. 
Therefore, BH-HMXBs in the early universe are an impor- 
tant ingredient in reconciling the apparent disparity be- 
tween the observed number of dwarf galaxies in the Galactic 
halo with the number of low-mass galaxies predicted by the 
cold dark matter model of the universe. 



9. Conclusions 

The main conclusions of this work are the following: 

1. The ratio of black holes to neutron stars and the ratio 
of black hole binaries to solitary black holes should in- 
crease with redshift; that is, the rate of formation of BH- 
HMXBs was significantly higher in the early Universe 
than at present. 

2. Feedback from one of those BH-HMXBs during its 
whole lifetime can be more than 10 54 erg (orders of mag- 
nitude larger than that of the photonic and barionic 
energy from a typical core collapse supernova). 

3. An accreting black hole in a high- mass binary emits a 
total number of ionizing photons that is comparable to 
its progenitor star, but one X-ray photon emitted by an 
accreting black hole may cause the ionization of several 
tens of hydrogen atoms in a fully neutral medium. 

4. The most important effect of BH-HMXBs in the early 
universe is the heating of the IGM. Soft X-rays and 
inverse-Compton scattering from relativistic electrons 
produced by BH-HMXBs heat the low-density medium 
over large distances to temperatures of ss 10 4 K, which 
limits the recombination rate of hydrogen and keeps the 
IGM ionized. 



5. A temperature of the IGM of ss 10 4 K limits the forma- 
tion of faint galaxies at high redshifts. It constrains the 
total mass of dwarf galaxies to > 10 9 Mq . 

6. BH-HMXBs in the early universe are important ingredi- 
ents for reconciling the apparent disparity between the 
observed number of faint dwarf galaxies in the Galactic 
halo with the number of low-mass galaxies predicted by 
the cold dark matter model of the universe. 

7. An additional effect of metallicity in the formation of 
BH-HMXBs (jMirabel 2010|) is to boost the formation of 
BH-BH binaries as more likely sources of gravitational 
waves than NS-NS systems ( |Belczynski et al. 2010) . 
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